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Abstract 

Astrophysical and cosmological observations do not require the dark matter parti- 
cles to be absolutely stable. If they are indeed unstable, their decay into positrons 
might occur at a sufficiently large rate to allow the indirect detection of dark mat- 
ter through an anomalous contribution to the cosmic positron flux. In this paper 
we discuss the implications of the excess in the positron fraction recently reported 
by the PAMELA collaboration for the scenario of decaying dark matter. To this 
end, we have performed a model-independent analysis of possible signatures by 
studying various decay channels in the case of both a fermionic and a scalar dark 
matter particle. We find that the steep rise in the positron fraction measured 
by PAMELA at energies larger than 10 GeV can naturally be accommodated in 
several realizations of the decaying dark matter scenario. 
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1 Introduction 



The spallation of primary cosmic-ray protons and other nuclei on the interstellar 
medium produces a flux of secondary positrons which is expected to decrease monoton- 
ically with the energy [1]. Interestingly, measurements of cosmic-ray positrons under- 
taken by a series of experiments over the last twenty years, HEAT [2], CAPRICE [3], 
MASS [4] and AMS-01 [5], indicated the existence of an excess of positrons at energies 
above 7 GeV with respect to the expectations from a purely secondary component. 
The discovery of this excess raised a lot of interest among the particle physics and 
astrophysics communities, which interpreted the excess as a possible indirect signature 
of dark matter. 

If dark matter particles are weakly interacting, they annihilate in the center of our 
Galaxy, producing a primary flux of positrons. This possibility has been extensively 
discussed over the last years as a potential explanation of the HEAT anomaly, although 
typically large boost factors have to be invoked in order to achieve sufficiently high 
fluxes [6-8]. 

Nevertheless, dark matter self-annihilation is not the only possibility for the in- 
direct detection of the dark matter. Strictly speaking, the viability of a particle as 
a dark matter candidate does not require its absolute stability, but merely that the 
dark matter lifetime is longer than the age of the Universe. Then, if the dark matter 
decays proceed at a sufficiently high rate, the decay products might be detectable. 
There are in fact some physically well-motivated dark matter candidates which decay 
with very long lifetimes. For instance, gravitino dark matter which is unstable due 
to a small breaking of i?-parity constitutes a very interesting scenario that leads to a 
thermal history of the Universe consistent with the observed abundances of the pri- 
mordial elements, the observed dark matter relic abundance and the observed baryon 
asymmetry [9]. The late decays of dark matter gravitinos produce a flux of gamma 
rays [10,11], positrons [12,11], antiprotons [12] and neutrinos [13] which contribute 
to the total fluxes received at the Earth. Remarkably, it has been pointed out that 
the EGRET anomaly in the extragalactic gamma-ray background [14] and the HEAT 
excess in the positron fraction [2] can be simultaneously explained by the decay of 
dark matter gravitinos with a mass m ~ 150 GeV and a lifetime r ~ 10^^ s [12,11]. 
Other candidates for decaying dark matter with electroweak masses include hidden 
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gauge bosons [15], hidden gauginos [16], right-handed sneutrinos in i?-parity breaking 
scenarios [17] or baryonic bound states of messenger quarks [18]. 

Very recently, the PAMELA collaboration [19] has pubhshed measurements of the 
cosmic-ray positron fraction performed with unprecedented accuracy [20] . These mea- 
surements have not only confirmed a significant deviation with respect to the expecta- 
tions from a purely secondary component, but have also provided evidence for a very 
sharp rise of the spectrum at energies 7 - 100 GeV. In view of the new results, it is 
important to study whether the scenario of decaying dark matter is consistent with the 
energy spectrum measured by PAMELA and what constraints the new data impose 
on the nature of decaying dark matter. It should be borne in mind, though, that the 
astrophysical uncertainties in the determination of the secondary positron component 
are still large [21] and that nearby astrophysical sources such as pulsars might produce 
sizable positron fluxes in the energy range explored by PAMELA [22]. 

In order to keep the analysis as model-independent as possible, we will analyze the 
cases that the dark matter particle is either a fermion or a scalar, and we will compute 
the predictions for the positron fraction for various decay channels and different dark 
matter masses and lifetimes. Namely, in the case of a fermionic dark matter particle 
ip, wc will consider the two-body decay channels ip Z^u.ip W^l^, as well as the 
three-body decay channels ip £~^£~i', with i = e, n, r being the charged leptons. 
On the other hand, for a scalar dark matter particle 0, we will consider the two-body 
decay channels, ^ Z^Z^, (f) W+W' , ^ 1+^. 

This paper is organized as follows: in Section 2 we will review the propagation 
of positrons in the Galaxy. In Section 3 we will present our results for the positron 
fraction expected from the decay of a fermionic or a scalar dark matter particle and 
we will discuss the sensitivity of the results to the choice of the propagation model. 
Lastly, in Section 4 we will present our conclusions. 

2 Positron Propagation 

Positron propagation in the Milky Way is commonly described by a stationary two-zone 
diffusion model with cylindrical boundary conditions [23]. Under this approximation, 
the number density of positrons per unit energy, fg+{E,f,t), satisfies the following 
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transport equation: 

= ^ = V • [K{E,r)VU] + -^[h{E,T)U] + Q^E.t) , (1) 

where convection and annihilations in the Galactic disk are neglected. The boundary 
conditions require the solution fg+ {E,f,t) to vanish at the boundary of the diffusion 
zone, which is approximated by a cylinder with half-height L — 1 — 15 kpc and radius 
i? = 20 kpc. 

The first term on the right-hand side of the transport equation is the diffusion 
term, which accounts for the propagation of positrons through the tangled Galactic 
magnetic fields. The diffusion coefficient K[E, r) is assumed to be constant throughout 
the diffusion zone and is parametrized by: 

K{E) ^KoPn' , (2) 

where (5 — v/c is the velocity and TZ is the rigidity of the particle, which is defined 
as the momentum in GeV per unit charge, TZ = p{GeV)/Z. The normalization Kq 
and the spectral index 5 of the diffusion coefficient are related to the properties of the 
interstellar medium and can be determined from flux measurements of other cosmic- 
ray species, mainly from the Boron-to-Garbon (B/G) ratio [24]. We list in Table 1 
the diffusion parameters for the propagation models M2, MED and Ml proposed in 
[25], which are consistent with the observed B/G ratio. The second term accounts for 
energy losses due to inverse Gompton scattering on starlight and the cosmic microwave 
background, as well as synchrotron radiation and ionization. The rate of energy loss, 
b{E,r), is assumed to be a spatially constant function parametrized by b{E) — 
with Eq — 1 GeV and te — 10^^ s. Lastly, Qe+{E,r) is the source term of positrons 
from the decay of a dark matter particle with mass ttidm and lifetime tdm' 

Q^,^E,f, = J^^-M^, (3) 

mBMTDM aE 

where dN^+/dE is the energy spectrum of positrons produced in the decay and pDM(r) 
is the density profile of dark matter in our Galaxy. For our numerical analysis, we will 
adopt the spherically symmetric Navarro, Frenk and White (NEW) profile [26]: 

PDM(r) = , (4) 

with po = 0.26 GeV/cm^ and Tc = 20 kpc. The normalization is chosen such that the 
local halo density is PdmI?^©) — 0.3 GeV/cm^ with = 8.5 kpc [27]. 
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Model 


S 


Ko (kpcVMyr) 


L (kpc) 


M2 


0.55 


0.00595 


1 


MED 


0.70 


0.0112 


4 


Ml 


0.46 


0.0765 


15 



Table 1: Diffusion parameters for the propagation models M2, MED and Ml proposed in 
[25] which are consistent with the observed B/C ratio. 



The solution of the transport equation at the Solar System, r = Tq, z = 0, can be 
formally expressed by the convolution 

feAE) = dE'GAE. E')^-^^^ . (5) 

"^DM^DM Jo ah' 

The solution is thus factorized into two parts. The first part, given by the Green's 
function G{E,E'), encodes all of the information about the astrophysics (such as the 
details of the halo profile and the complicated propagation of positrons in the Galaxy) 
and is universal for any decaying dark matter candidate. The remaining part depends 
exclusively on the nature and properties of the decaying dark matter candidate, namely 
the mass, the lifetime and the energy spectrum of positrons produced in the decay. In 
the next Section we will analyze several phenomenological scenarios of decaying dark 
matter and discuss their viability in view of the PAMELA data. 
The explicit form of the Green's function is [8] 
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Ge^E, E') = J2 Bnm(E, E')Jo (Cn^) sm (^) , (6) 

n,m=l 

where Jq is the zeroth-order Bessel function of the first kind, whose successive zeros 
are denoted by C„. On the other hand. 



Bnm{E,E') = ^l^CnrnGXp 

with 




eV'^ fE'^^-^ 



EqJ \Eq 



(7) 



2 f f f r'\ [rmr i 

^ iKcnm Jo L [^-R ) [ 2L - "'^J ' 

where Ji is the first-order Bessel function. 

The Green's function can be well approximated by the following interpolating func- 
tion, which is valid for any decaying dark matter particle [12] : 

1 

G.^E, E') ~ :^e"+^(^'"'-^"")^(E' - E) cm'^ s , (9) 



model 


a b 


M2 
MED 
Ml 


-0.9716 -10.012 

-1.0203 -1.4493 
-0.9809 -1.1456 



Table 2: Coefficients of the interpolating function Eq. (9) for the positron Green's function, 
assuming a NFW halo profile and for the different diffusion models in Table 1. 

where E and E' arc expressed in units of GeV. The coefficients a and b can be found 
in Table 2 for the NFW profile and the different diffusion models listed in Table 1. 
This approximation works better than 15 - 20% over the whole range of energies. We 
find numerically that the Green's function is not very sensitive to the choice of the 
halo profile, since the Earth receives only positrons created within a few kpc from the 
Sun, where the different halo profiles are very similar. Therefore, the Green's function 
coefficients for other halo profiles can be well approximated by the values in Table 2. 

Finally, the fiux of primary positrons at the Solar System from dark matter decay 
is given by: 

^TiE) = ^fem- (10) 

In order to compare our results with experiments, we will calculate the positron 
fraction, which is defined as the ratio of the total positron fiux over the total electron 
plus positron fiuxes, $e+/(^e- + ^e+)- The total positron fiux receives contributions 
from the dark matter decay as well as from a secondary component stemming from the 
collision of primary protons and other nuclei on the interstellar medium, which consti- 
tutes the background to any dark matter signal. On the other hand, the total electron 
flux has a primary astrophysical component, presumably originating from supernova 
remnants, a secondary component from spallation of cosmic rays on the interstellar 
medium, and an exotic primary component from dark matter decay that might be im- 
portant at high energies in some scenarios. For the background fiuxes of primary and 
secondary electrons, as well as secondary positrons, we will use the parametrizations 
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obtained in [7] from detailed computer simulations of cosmic-ray propagation [1]: 

= ^^^^ — (GeV-^cm-'s-^sr-i) , (11) 

" ^ ' 1 + 11 + 3.2^2.15 ^ ^ ' y > 

n 70 F^-'^ 

*~<^' = l + 110^.-^ + 600^^» + 580E^.^ (GeV-W-Vsr-) . (12) 

^fHE) = ^^:!-=^ — (GeV-^cm-V^sr-i) , (13) 

' 1 + 650 + 1500 £;4-2 ^ ' ' ^ ' 

where E is the energy expressed in units of GeV. 



3 Results 

For a given model of decaying dark matter and a given propagation model, defined 
by the parameters in Table 1, the predictions for the positron fraction at Earth can 
be readily computed using the formalism explained in the previous section. To derive 
our results wc solved the transport equation Eq. (1) numerically, although for a quick 
estimate of the positron fraction the interpolating function Eq. (9) may also be used. 

To keep the analysis as model-independent as possible, we will analyze several 
possibilities for the decaying dark matter, computing the prediction for the positron 
fraction for either a fermionic or a bosonic particle which decays in various channels 
with a branching ratio of 100%. We will treat the dark matter mass and lifetime as 
free parameters, while the normalization of the background is kept fixed. We will first 
adopt the MED propagation model in order to better compare the predictions from 
different particle physics scenarios and later on, we will analyze the sensitivity of the 
predictions to the choice of the propagation model. 

Let us now discuss the cases of fermionic and scalar dark matter particles separately. 



In the case that the dark matter particle is a fermion -0, the following decay channels 
are possible: 

i) £+ru , (14) 

provided that the respective channels arc kincmatically open. 

The fragmentation of the boson in the decay Z^v produces a continu- 

ous spectrum of positrons (mainly from tt"*" decay) that we have obtained using the 
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Figure 1 : Positron fraction from the decay of the fermionic dark matter particle in the channel 
tp — > Z^u when the dark matter mass is, from left to right, moM = 150, 300, 600, 1000 GeV. 
The lifetime has been chosen to provide a qualitatively good fit to the data and is ~ 5 x 1025 s 
in all the cases. 

event generator PYTHIA 6.4 [28]^. The predicted positron fraction is shown in Fig. 1, 
compared to the PAMELA and HEAT data, for the MED propagation model and 
for different dark matter masses. The lifetime has been chosen in order to produce a 
qualitatively good fit of the prediction to the data points. We operate under the as- 
sumption that the difference between the low-energy data from HEAT and PAMELA 
is due to solar modulation. Therefore, we use the HEAT data at low energies, which 
were recorded during a period of minimum solar activity. It is apparent from the figure 
that this decay channel by itself cannot explain the steep rise of the spectrum observed 
by PAMELA for any of the masses analyzed here. We have also checked that using 
different propagation models does not improve the fit to the data significantly. 

On the other hand, we show in Fig. 2 the prediction for the positron fraction 
when the fermionic dark matter particle decays as ^ Wi. The positrons created 
in the fragmentation of the W gauge bosons produce a rather fiat contribution to the 
positron fraction. However, the hard positrons resulting from the decay of the fi and 
T leptons or directly from the dark matter decay into positrons produce a rise in the 
spectrum, which is most prominent in the decay mode ip — > W^e^, although it is also 

"'^The fragmentation of the weak gauge bosons also produces fluxes of primary antiprotons and 
gamma-rays, which are severely constrained by the PAMELA [29] and EGRET [30] experiments, 
respectively. The predictions for the antiproton and the gamma-ray fluxes of the phenomenological 
scenarios studied in this paper wiU be presented elsewhere [31]. 
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Figure 2: Positron fraction from the decay of the fermionic dark matter particle in the 
channels ip W^e^ (top-left panel), W^/j,^ (top-right panel) and ip —>■ W^t^ (bottom 
panel), when the dark matter mass is, from left to right, mDM = 150, 300, 600, 1000 GeV. 
The lifetime, which ranges between lO'^^ s and 5 x 10^^ s, is different in each case and has 
been chosen to provide a qualitatively good fit to the data. 

quite visible in ^ W^fi^. Interestingly, these two decay models can qualitatively 
reproduce the energy spectrum measured by PAMELA for dark matter masses larger 
than ~ 300 GeV. Future measurements of the positron fraction at energies 100 - 300 
GeV, as planned by the PAMELA collaboration, will be crucial to discriminate among 
these possibilities. Note that when the dark matter mass is large, the positron fraction 
in the scenario with direct decay into positrons, ip W^e^, can reach values as large 
as 30 - 40%. In this case, not only the positron flux, but also the electron flux will 
receive a significant primary contribution from dark matter decay. 

Lastly, we have also analyzed the case that the dark matter particle decays lepton- 
ically in a three-body decay, tp —> i'^i'u; the results are shown in Fig. 3. The decay 
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Figure 3: Positron fraction from the decay of the fermionic dark matter particle in the 
channels ^ — > e^e~z^ (top-left panel), ip — > fi~^fi~i' (top-right panel) and ip — > t~^t~i' (bottom 
panel), when the dark matter mass is, from left to right, ttt-dm = 150, 300, 600, 1000 GeV. 
The lifetime, which ranges between 5 x 10^^ s and 8 x 10^^ s, is different in each case and has 
been chosen to provide a qualitatively good fit to the data. 

-0 e^e^v can explain the PAMELA anomaly when the dark matter mass is larger 
than ~ 300 GeV, while the decay channel if) —> requires larger masses. In 

contrast, the energy spectrum produced in the decay ip t^t~i>, although it exhibits 
a notable bump, is too flat to explain the observations by itself. In this case, other 
contributions to the positron flux should be invoked, for instance from pulsars, in order 
to reproduce the observations. 



When the decaying dark matter particle is a scalar, the following decay modes are 
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Figure 4: Positron fraction from the decay of a scalar dark matter particle in the channels 
(j) — > Z^Z^ (top-left panel) and (j) W~^W~ (top-right panel) when the dark matter mass 
is, from left to right, m-DM = 300, 600, 1000 GeV. The lifetime has been chosen to provide 
a qualitatively good fit to the data and is ~ 2 x 10^^ s (~ 10^^ s) for the decay into Z {W) 
bosons. 

possible: 

(f) i+r, (15) 

again provided that the decays are kinematically open. 

We show in Fig. 4 the expected positron fraction from the decay of a scalar particle 
into weak gauge bosons, Z^Z^ or — W^W~ . In these decay modes, no hard 
lepton is produced. Instead, only positrons from the fragmentation of the weak gauge 
bosons will contribute. As a result, the spectral shape of the positron fraction is too 
flat to explain the steep rise observed by PAMELA by itself. 

We have also calculated the predictions for the positron fraction when the scalar 
dark matter decays directly into two charged leptons, £~^£~; the results are shown 
in Fig. 5. The predictions are qualitatively very similar to the ones in the scenario 
where a dark matter fermion decays into W^i"^, the most noticeable difference be- 
ing the larger cutoff in the energy spectrum in the scalar case, for the same dark 
matter mass. To be more precise, in the scenario with ip W^i"^ the cutoff is at 
^^01^1/2(1 — M^/m^jy[) while in the scenario with it is at the larger value 
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Figure 5: Positron fraction from the decay of a scalar dark matter particle in the channels 
(j) — > e"'"e~ (top-left panel), 4) — > ^'^ ii~ (top-right panel) and (j) — > r^r~ (bottom panel), 
when the dark matter mass is, from left to right, mDM = 150, 300, 600, 1000 GeV. The 
lifetime, which ranges between 10^^ s and 10^^ s, is different in each case and has been chosen 
to provide a qualitatively good fit to the data. 

For this class of scenarios, the decay modes into positrons or antimuons can 
accommodate the PAMELA anomaly better than the decay mode into antitaus. 



Above we have analyzed the predictions for the positron fraction for different sce- 
narios of fermionic or scalar decaying dark matter, assuming the MED propagation 
model. After discussing the sensitivity of the predictions to the choice of the par- 
ticle physics model, we would like to briefly address the sensitivity of the results to 
astrophysical uncertainties. 

We show in Fig. 6 the sensitivity to the propagation model of four particle physics 
scenarios consisting of a dark matter fermion which decays into two particles, ip —>■ 
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W'^le, iJ,)^, or into three particles, ip (e, //)+(e, These scenarios are charac- 
terized by the direct decay into hard positrons or antimuons, and were found to be 
compatible with the excess observed by PAMELA for the MED propagation model. 
The results in Fig. 6 show that this conclusion still holds for the Ml and M2 propa- 
gation models: the Ml and MED propagation models yield very similar predictions of 
the positron fraction, while the M2 propagation model yields a slightly larger positron 
fraction at high energies and a smaller positron fraction at low energies. The rea- 
son is that in the M2 propagation model the diffusion is minimal and thus the hard 
spectrum of injected positrons from dark matter decay is less altered by the propa- 
gation than in the MED and Ml model, yielding a steeper rise in the positron frac- 
tion measured at Earth. The same conclusion holds for the scenarios where a scalar 
dark matter particle directly decays into an electron-positron pair or a muon-antimuon 
pair, (f) {e, ii)'^{e, , which as discussed above yield very similar signatures in the 
positron fraction as the scenarios with fermionic dark matter ip — > W^le, /i)^, where 
the dark matter also decays into a monoenergetic positron or antimuon. 

As a last remark, let us note that the origin of the positron excess might not be the 
decay of the dark matter itself, but the decay of another long-lived particle which is 
present in our Galaxy. Notice that in the scenarios analyzed in this paper with decaying 
dark matter, the shape of the spectrum is determined by the dark matter mass, while 
the normalization, by the combination pDMl^'oj/TDM, with pdm(^0) = 0.3 GeV/cm^ 
being the local density of dark matter particles (c/. Eqs. (3,4)). Therefore, identical 
signatures in the positron fraction are obtained if the decaying particle has a mass m', 
a hfetime r' and a local abundance p'{rQ) satisfying^ 

m' = moM , 

p'{rQ)/T' = PDM(r0)/TDM ■ (16) 

This opens new possibihties for model building. For instance, it is conceivable that the 
dominant component of dark matter is completely stable (perhaps even annihilating), 
while the origin of the positron excess is the decay of a subdominant form of dark 
matter. 

^We assume here that any weakly or superweakly interacting particle with a lifetime longer than 
the age of the Universe is present in our Galaxy and has a halo distribution identical to the dominant 
component of dark matter. 
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Figure 6: Dependence of the predicted positron fraction on the propagation model. For a 
decaying dark matter fermion with ttt-dm = 1000 GeV, we show the decay modes ^ —>■ W^e^ 
(top-left panel), ip W^fj,^ (top-right panel) ijj — > e'^e'v (bottom-left panel), Tp fi'^fi^v 
(bottom-right panel). The dotted line corresponds to the M2 propagation model, while the 
almost indistinguishable dashed-dotted lines correspond to the MED and the Ml models (see 
Table 1). 

4 Conclusions 

We have computed the predictions for the positron fraction of several decaying dark 
matter scenarios and discussed their viability in view of the new measurements reported 
by the PAMELA collaboration. We have studied scenarios where the dark matter 
particle is a fermion ip, which decays in the two-body decay channels ip Z^u, W^i"^ 
or in the three-body decay channel ip z/, where i = e, n, t denotes the charged 

leptons. In addition, we have also studied scenarios where the dark matter particle is 
a scalar (p, which decays as ^ (f) W'^W~ , (p — * 
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We have found that decay channels producing hard positrons and antimuons are 
favored by the PAMELA data. Namely, the decay channels of a fermionic dark matter 
particle — > W{ii,e) and ip — > e)'^(/x, e)~i/, and the decay channels of a scalar 
dark matter particle — > (//, e) e)~ produce a steep rise in the positron fraction 
above 10 GeV, as observed by the PAMELA collaboration. It is worth noting that 
the spectral shape of the positron fraction resulting from a sharply peaked positron 
or antimuon injection spectrum matches the experimental results quite accurately. On 
the other hand, decay modes involving the tau lepton can contribute to an excess in the 
positron fraction at high energies, although the resulting spectrum is too flat. Finally, 
the decay of a scalar particle into weak gauge bosons, (j) Z^Z^, 4> — > W~^W~, can 
also contribute to an excess in the positron fraction, although the spectrum is also 
too flat to fit the PAMELA data by itself. In addition, these decay channels are also 
disfavored by antiproton overproduction constraints for dark matter lifetimes that yield 
appreciable contributions to the positron flux. 

If dark matter decay is indeed the primary cause for the observed positron excess, 
the PAMELA data clearly seem to point toward a rather heavy dark matter particle, 
^DM ^ 300 GeV, which preferentially decays directly into first or second generation 
charged Icptons with a lifetime tdm ~ 10^^ s. 

Future measurements of the positron fraction in the energy range 100 - 300 GeV 
by PAMELA will provide invaluable information about scenarios with decaying dark 
matter and the properties of the decaying particles. Furthermore, the constraints on 
the predicted antiproton flux by these same scenarios from PAMELA [29] and on the 
gamma-ray flux from EGRET [30] (and in the near future from the Fermi Gamma-ray 
Space Telescope [32]) offer complementary information about this scenario, which will 
be presented elsewhere [31]. 

Note Added 

During the completion of this work two preprints appeared presenting related analyses 
[33,34]. 
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